home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_tbsv.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  54.9 KB  |  1,820 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_tbsv (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int order = 101;
  13.    int trans = 111;
  14.    int uplo = 121;
  15.    int diag = 131;
  16.    int N = 3;
  17.    int K = 1;
  18.    int lda = 3;
  19.    float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
  20.    float X[] = { -0.305f, 0.61f, -0.831f };
  21.    int incX = -1;
  22.    float x_expected[] = { -0.354651f, -2.40855f, 0.481076f };
  23.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  24.    {
  25.      int i;
  26.      for (i = 0; i < 3; i++) {
  27.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1230)");
  28.      }
  29.    };
  30.   };
  31.  
  32.  
  33.   {
  34.    int order = 101;
  35.    int trans = 111;
  36.    int uplo = 121;
  37.    int diag = 132;
  38.    int N = 3;
  39.    int K = 1;
  40.    int lda = 3;
  41.    float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
  42.    float X[] = { -0.305f, 0.61f, -0.831f };
  43.    int incX = -1;
  44.    float x_expected[] = { -0.305f, 0.84973f, -1.00859f };
  45.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  46.    {
  47.      int i;
  48.      for (i = 0; i < 3; i++) {
  49.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1231)");
  50.      }
  51.    };
  52.   };
  53.  
  54.  
  55.   {
  56.    int order = 101;
  57.    int trans = 111;
  58.    int uplo = 122;
  59.    int diag = 131;
  60.    int N = 3;
  61.    int K = 1;
  62.    int lda = 3;
  63.    float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
  64.    float X[] = { -0.305f, 0.61f, -0.831f };
  65.    int incX = -1;
  66.    float x_expected[] = { -2.71619f, -1.09055f, -3.97608f };
  67.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  68.    {
  69.      int i;
  70.      for (i = 0; i < 3; i++) {
  71.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1232)");
  72.      }
  73.    };
  74.   };
  75.  
  76.  
  77.   {
  78.    int order = 101;
  79.    int trans = 111;
  80.    int uplo = 122;
  81.    int diag = 132;
  82.    int N = 3;
  83.    int K = 1;
  84.    int lda = 3;
  85.    float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
  86.    float X[] = { -0.305f, 0.61f, -0.831f };
  87.    int incX = -1;
  88.    float x_expected[] = { -0.56589f, 0.303361f, -0.831f };
  89.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  90.    {
  91.      int i;
  92.      for (i = 0; i < 3; i++) {
  93.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1233)");
  94.      }
  95.    };
  96.   };
  97.  
  98.  
  99.   {
  100.    int order = 102;
  101.    int trans = 111;
  102.    int uplo = 121;
  103.    int diag = 131;
  104.    int N = 3;
  105.    int K = 1;
  106.    int lda = 3;
  107.    float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
  108.    float X[] = { -0.305f, 0.61f, -0.831f };
  109.    int incX = -1;
  110.    float x_expected[] = { 1.30901f, -0.656172f, -5.13458f };
  111.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  112.    {
  113.      int i;
  114.      for (i = 0; i < 3; i++) {
  115.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1234)");
  116.      }
  117.    };
  118.   };
  119.  
  120.  
  121.   {
  122.    int order = 102;
  123.    int trans = 111;
  124.    int uplo = 121;
  125.    int diag = 132;
  126.    int N = 3;
  127.    int K = 1;
  128.    int lda = 3;
  129.    float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
  130.    float X[] = { -0.305f, 0.61f, -0.831f };
  131.    int incX = -1;
  132.    float x_expected[] = { -0.305f, 0.8723f, -0.509121f };
  133.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  134.    {
  135.      int i;
  136.      for (i = 0; i < 3; i++) {
  137.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1235)");
  138.      }
  139.    };
  140.   };
  141.  
  142.  
  143.   {
  144.    int order = 102;
  145.    int trans = 111;
  146.    int uplo = 122;
  147.    int diag = 131;
  148.    int N = 3;
  149.    int K = 1;
  150.    int lda = 3;
  151.    float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
  152.    float X[] = { -0.305f, 0.61f, -0.831f };
  153.    int incX = -1;
  154.    float x_expected[] = { 0.524539f, -0.961964f, 1.22026f };
  155.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  156.    {
  157.      int i;
  158.      for (i = 0; i < 3; i++) {
  159.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1236)");
  160.      }
  161.    };
  162.   };
  163.  
  164.  
  165.   {
  166.    int order = 102;
  167.    int trans = 111;
  168.    int uplo = 122;
  169.    int diag = 132;
  170.    int N = 3;
  171.    int K = 1;
  172.    int lda = 3;
  173.    float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
  174.    float X[] = { -0.305f, 0.61f, -0.831f };
  175.    int incX = -1;
  176.    float x_expected[] = { -0.920972f, 0.783679f, -0.831f };
  177.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  178.    {
  179.      int i;
  180.      for (i = 0; i < 3; i++) {
  181.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1237)");
  182.      }
  183.    };
  184.   };
  185.  
  186.  
  187.   {
  188.    int order = 101;
  189.    int trans = 112;
  190.    int uplo = 121;
  191.    int diag = 131;
  192.    int N = 3;
  193.    int K = 1;
  194.    int lda = 3;
  195.    float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
  196.    float X[] = { 0.144f, 0.635f, 0.116f };
  197.    int incX = -1;
  198.    float x_expected[] = { 16.8676f, 17.3503f, 5.27273f };
  199.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  200.    {
  201.      int i;
  202.      for (i = 0; i < 3; i++) {
  203.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1238)");
  204.      }
  205.    };
  206.   };
  207.  
  208.  
  209.   {
  210.    int order = 101;
  211.    int trans = 112;
  212.    int uplo = 121;
  213.    int diag = 132;
  214.    int N = 3;
  215.    int K = 1;
  216.    int lda = 3;
  217.    float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
  218.    float X[] = { 0.144f, 0.635f, 0.116f };
  219.    int incX = -1;
  220.    float x_expected[] = { 0.209676f, 0.54278f, 0.116f };
  221.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  222.    {
  223.      int i;
  224.      for (i = 0; i < 3; i++) {
  225.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1239)");
  226.      }
  227.    };
  228.   };
  229.  
  230.  
  231.   {
  232.    int order = 101;
  233.    int trans = 112;
  234.    int uplo = 122;
  235.    int diag = 131;
  236.    int N = 3;
  237.    int K = 1;
  238.    int lda = 3;
  239.    float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
  240.    float X[] = { 0.144f, 0.635f, 0.116f };
  241.    int incX = -1;
  242.    float x_expected[] = { 0.212077f, -5.01482f, -1.14722f };
  243.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  244.    {
  245.      int i;
  246.      for (i = 0; i < 3; i++) {
  247.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1240)");
  248.      }
  249.    };
  250.   };
  251.  
  252.  
  253.   {
  254.    int order = 101;
  255.    int trans = 112;
  256.    int uplo = 122;
  257.    int diag = 132;
  258.    int N = 3;
  259.    int K = 1;
  260.    int lda = 3;
  261.    float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
  262.    float X[] = { 0.144f, 0.635f, 0.116f };
  263.    int incX = -1;
  264.    float x_expected[] = { 0.144f, 0.615848f, 0.242249f };
  265.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  266.    {
  267.      int i;
  268.      for (i = 0; i < 3; i++) {
  269.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1241)");
  270.      }
  271.    };
  272.   };
  273.  
  274.  
  275.   {
  276.    int order = 102;
  277.    int trans = 112;
  278.    int uplo = 121;
  279.    int diag = 131;
  280.    int N = 3;
  281.    int K = 1;
  282.    int lda = 3;
  283.    float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
  284.    float X[] = { 0.144f, 0.635f, 0.116f };
  285.    int incX = -1;
  286.    float x_expected[] = { 1.28844f, -5.49514f, 0.145912f };
  287.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  288.    {
  289.      int i;
  290.      for (i = 0; i < 3; i++) {
  291.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1242)");
  292.      }
  293.    };
  294.   };
  295.  
  296.  
  297.   {
  298.    int order = 102;
  299.    int trans = 112;
  300.    int uplo = 121;
  301.    int diag = 132;
  302.    int N = 3;
  303.    int K = 1;
  304.    int lda = 3;
  305.    float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
  306.    float X[] = { 0.144f, 0.635f, 0.116f };
  307.    int incX = -1;
  308.    float x_expected[] = { 0.0563823f, 0.65878f, 0.116f };
  309.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  310.    {
  311.      int i;
  312.      for (i = 0; i < 3; i++) {
  313.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1243)");
  314.      }
  315.    };
  316.   };
  317.  
  318.  
  319.   {
  320.    int order = 102;
  321.    int trans = 112;
  322.    int uplo = 122;
  323.    int diag = 131;
  324.    int N = 3;
  325.    int K = 1;
  326.    int lda = 3;
  327.    float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
  328.    float X[] = { 0.144f, 0.635f, 0.116f };
  329.    int incX = -1;
  330.    float x_expected[] = { 1.08271f, -3.73662f, 140.301f };
  331.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  332.    {
  333.      int i;
  334.      for (i = 0; i < 3; i++) {
  335.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1244)");
  336.      }
  337.    };
  338.   };
  339.  
  340.  
  341.   {
  342.    int order = 102;
  343.    int trans = 112;
  344.    int uplo = 122;
  345.    int diag = 132;
  346.    int N = 3;
  347.    int K = 1;
  348.    int lda = 3;
  349.    float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
  350.    float X[] = { 0.144f, 0.635f, 0.116f };
  351.    int incX = -1;
  352.    float x_expected[] = { 0.144f, 0.652424f, -0.402677f };
  353.    cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  354.    {
  355.      int i;
  356.      for (i = 0; i < 3; i++) {
  357.        gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1245)");
  358.      }
  359.    };
  360.   };
  361.  
  362.  
  363.   {
  364.    int order = 101;
  365.    int trans = 111;
  366.    int uplo = 121;
  367.    int diag = 131;
  368.    int N = 3;
  369.    int K = 1;
  370.    int lda = 3;
  371.    double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
  372.    double X[] = { 0.332, 0.588, 0.252 };
  373.    int incX = -1;
  374.    double x_expected[] = { -0.967930029155, 0.138412575592, 0.506166027443 };
  375.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  376.    {
  377.      int i;
  378.      for (i = 0; i < 3; i++) {
  379.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1246)");
  380.      }
  381.    };
  382.   };
  383.  
  384.  
  385.   {
  386.    int order = 101;
  387.    int trans = 111;
  388.    int uplo = 121;
  389.    int diag = 132;
  390.    int N = 3;
  391.    int K = 1;
  392.    int lda = 3;
  393.    double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
  394.    double X[] = { 0.332, 0.588, 0.252 };
  395.    int incX = -1;
  396.    double x_expected[] = { 0.332, 0.819736, 0.615143048 };
  397.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  398.    {
  399.      int i;
  400.      for (i = 0; i < 3; i++) {
  401.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1247)");
  402.      }
  403.    };
  404.   };
  405.  
  406.  
  407.   {
  408.    int order = 101;
  409.    int trans = 111;
  410.    int uplo = 122;
  411.    int diag = 131;
  412.    int N = 3;
  413.    int K = 1;
  414.    int lda = 3;
  415.    double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
  416.    double X[] = { 0.332, 0.588, 0.252 };
  417.    int incX = -1;
  418.    double x_expected[] = { -0.364842154056, -0.326531140246, -0.568848758465 };
  419.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  420.    {
  421.      int i;
  422.      for (i = 0; i < 3; i++) {
  423.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1248)");
  424.      }
  425.    };
  426.   };
  427.  
  428.  
  429.   {
  430.    int order = 101;
  431.    int trans = 111;
  432.    int uplo = 122;
  433.    int diag = 132;
  434.    int N = 3;
  435.    int K = 1;
  436.    int lda = 3;
  437.    double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
  438.    double X[] = { 0.332, 0.588, 0.252 };
  439.    int incX = -1;
  440.    double x_expected[] = { 0.588397988, 0.747516, 0.252 };
  441.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  442.    {
  443.      int i;
  444.      for (i = 0; i < 3; i++) {
  445.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1249)");
  446.      }
  447.    };
  448.   };
  449.  
  450.  
  451.   {
  452.    int order = 102;
  453.    int trans = 111;
  454.    int uplo = 121;
  455.    int diag = 131;
  456.    int N = 3;
  457.    int K = 1;
  458.    int lda = 3;
  459.    double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
  460.    double X[] = { 0.332, 0.588, 0.252 };
  461.    int incX = -1;
  462.    double x_expected[] = { -0.550580431177, -0.571849444278, 0.248263427151 };
  463.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  464.    {
  465.      int i;
  466.      for (i = 0; i < 3; i++) {
  467.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1250)");
  468.      }
  469.    };
  470.   };
  471.  
  472.  
  473.   {
  474.    int order = 102;
  475.    int trans = 111;
  476.    int uplo = 121;
  477.    int diag = 132;
  478.    int N = 3;
  479.    int K = 1;
  480.    int lda = 3;
  481.    double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
  482.    double X[] = { 0.332, 0.588, 0.252 };
  483.    int incX = -1;
  484.    double x_expected[] = { 0.332, 0.701876, 0.696287508 };
  485.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  486.    {
  487.      int i;
  488.      for (i = 0; i < 3; i++) {
  489.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1251)");
  490.      }
  491.    };
  492.   };
  493.  
  494.  
  495.   {
  496.    int order = 102;
  497.    int trans = 111;
  498.    int uplo = 122;
  499.    int diag = 131;
  500.    int N = 3;
  501.    int K = 1;
  502.    int lda = 3;
  503.    double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
  504.    double X[] = { 0.332, 0.588, 0.252 };
  505.    int incX = -1;
  506.    double x_expected[] = { 1.50217883761, -1.21382140588, 0.407108239095 };
  507.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  508.    {
  509.      int i;
  510.      for (i = 0; i < 3; i++) {
  511.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1252)");
  512.      }
  513.    };
  514.   };
  515.  
  516.  
  517.   {
  518.    int order = 102;
  519.    int trans = 111;
  520.    int uplo = 122;
  521.    int diag = 132;
  522.    int N = 3;
  523.    int K = 1;
  524.    int lda = 3;
  525.    double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
  526.    double X[] = { 0.332, 0.588, 0.252 };
  527.    int incX = -1;
  528.    double x_expected[] = { 0.820345928, 0.699636, 0.252 };
  529.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  530.    {
  531.      int i;
  532.      for (i = 0; i < 3; i++) {
  533.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1253)");
  534.      }
  535.    };
  536.   };
  537.  
  538.  
  539.   {
  540.    int order = 101;
  541.    int trans = 112;
  542.    int uplo = 121;
  543.    int diag = 131;
  544.    int N = 3;
  545.    int K = 1;
  546.    int lda = 3;
  547.    double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
  548.    double X[] = { 0.58, 0.558, -0.54 };
  549.    int incX = -1;
  550.    double x_expected[] = { 18.994209959, 20.323927329, 2.7135678392 };
  551.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  552.    {
  553.      int i;
  554.      for (i = 0; i < 3; i++) {
  555.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1254)");
  556.      }
  557.    };
  558.   };
  559.  
  560.  
  561.   {
  562.    int order = 101;
  563.    int trans = 112;
  564.    int uplo = 121;
  565.    int diag = 132;
  566.    int N = 3;
  567.    int K = 1;
  568.    int lda = 3;
  569.    double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
  570.    double X[] = { 0.58, 0.558, -0.54 };
  571.    int incX = -1;
  572.    double x_expected[] = { 1.06925836, 0.72162, -0.54 };
  573.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  574.    {
  575.      int i;
  576.      for (i = 0; i < 3; i++) {
  577.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1255)");
  578.      }
  579.    };
  580.   };
  581.  
  582.  
  583.   {
  584.    int order = 101;
  585.    int trans = 112;
  586.    int uplo = 122;
  587.    int diag = 131;
  588.    int N = 3;
  589.    int K = 1;
  590.    int lda = 3;
  591.    double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
  592.    double X[] = { 0.58, 0.558, -0.54 };
  593.    int incX = -1;
  594.    double x_expected[] = { -3.27683615819, -4.47682615869, -1.97425326753 };
  595.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  596.    {
  597.      int i;
  598.      for (i = 0; i < 3; i++) {
  599.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1256)");
  600.      }
  601.    };
  602.   };
  603.  
  604.  
  605.   {
  606.    int order = 101;
  607.    int trans = 112;
  608.    int uplo = 122;
  609.    int diag = 132;
  610.    int N = 3;
  611.    int K = 1;
  612.    int lda = 3;
  613.    double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
  614.    double X[] = { 0.58, 0.558, -0.54 };
  615.    int incX = -1;
  616.    double x_expected[] = { 0.58, 0.11952, -0.53844624 };
  617.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  618.    {
  619.      int i;
  620.      for (i = 0; i < 3; i++) {
  621.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1257)");
  622.      }
  623.    };
  624.   };
  625.  
  626.  
  627.   {
  628.    int order = 102;
  629.    int trans = 112;
  630.    int uplo = 121;
  631.    int diag = 131;
  632.    int N = 3;
  633.    int K = 1;
  634.    int lda = 3;
  635.    double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
  636.    double X[] = { 0.58, 0.558, -0.54 };
  637.    int incX = -1;
  638.    double x_expected[] = { -6.6461072986, -0.788837290809, -1.78217821782 };
  639.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  640.    {
  641.      int i;
  642.      for (i = 0; i < 3; i++) {
  643.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1258)");
  644.      }
  645.    };
  646.   };
  647.  
  648.  
  649.   {
  650.    int order = 102;
  651.    int trans = 112;
  652.    int uplo = 121;
  653.    int diag = 132;
  654.    int N = 3;
  655.    int K = 1;
  656.    int lda = 3;
  657.    double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
  658.    double X[] = { 0.58, 0.558, -0.54 };
  659.    int incX = -1;
  660.    double x_expected[] = { 0.16345912, 0.55098, -0.54 };
  661.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  662.    {
  663.      int i;
  664.      for (i = 0; i < 3; i++) {
  665.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1259)");
  666.      }
  667.    };
  668.   };
  669.  
  670.  
  671.   {
  672.    int order = 102;
  673.    int trans = 112;
  674.    int uplo = 122;
  675.    int diag = 131;
  676.    int N = 3;
  677.    int K = 1;
  678.    int lda = 3;
  679.    double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
  680.    double X[] = { 0.58, 0.558, -0.54 };
  681.    int incX = -1;
  682.    double x_expected[] = { 0.767195767196, -82.9352869353, -123.564783625 };
  683.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  684.    {
  685.      int i;
  686.      for (i = 0; i < 3; i++) {
  687.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1260)");
  688.      }
  689.    };
  690.   };
  691.  
  692.  
  693.   {
  694.    int order = 102;
  695.    int trans = 112;
  696.    int uplo = 122;
  697.    int diag = 132;
  698.    int N = 3;
  699.    int K = 1;
  700.    int lda = 3;
  701.    double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
  702.    double X[] = { 0.58, 0.558, -0.54 };
  703.    int incX = -1;
  704.    double x_expected[] = { 0.58, 0.95124, -0.82822572 };
  705.    cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  706.    {
  707.      int i;
  708.      for (i = 0; i < 3; i++) {
  709.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1261)");
  710.      }
  711.    };
  712.   };
  713.  
  714.  
  715.   {
  716.    int order = 101;
  717.    int trans = 111;
  718.    int uplo = 121;
  719.    int diag = 131;
  720.    int N = 3;
  721.    int K = 1;
  722.    int lda = 3;
  723.    float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
  724.    float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
  725.    int incX = -1;
  726.    float x_expected[] = { 1.28871f, 0.289887f, 1.76043f, 1.27481f, 1.56506f, -2.35181f };
  727.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  728.    {
  729.      int i;
  730.      for (i = 0; i < 3; i++) {
  731.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1262) real");
  732.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1262) imag");
  733.      };
  734.    };
  735.   };
  736.  
  737.  
  738.   {
  739.    int order = 101;
  740.    int trans = 111;
  741.    int uplo = 121;
  742.    int diag = 132;
  743.    int N = 3;
  744.    int K = 1;
  745.    int lda = 3;
  746.    float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
  747.    float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
  748.    int incX = -1;
  749.    float x_expected[] = { 0.11f, 0.787f, -1.04259f, 0.18935f, 0.228474f, -0.564917f };
  750.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  751.    {
  752.      int i;
  753.      for (i = 0; i < 3; i++) {
  754.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1263) real");
  755.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1263) imag");
  756.      };
  757.    };
  758.   };
  759.  
  760.  
  761.   {
  762.    int order = 101;
  763.    int trans = 111;
  764.    int uplo = 122;
  765.    int diag = 131;
  766.    int N = 3;
  767.    int K = 1;
  768.    int lda = 3;
  769.    float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
  770.    float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
  771.    int incX = -1;
  772.    float x_expected[] = { -0.0906249f, 3.09442f, -1.60036f, 1.28475f, -0.582941f, 0.0383898f };
  773.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  774.    {
  775.      int i;
  776.      for (i = 0; i < 3; i++) {
  777.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1264) real");
  778.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1264) imag");
  779.      };
  780.    };
  781.   };
  782.  
  783.  
  784.   {
  785.    int order = 101;
  786.    int trans = 111;
  787.    int uplo = 122;
  788.    int diag = 132;
  789.    int N = 3;
  790.    int K = 1;
  791.    int lda = 3;
  792.    float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
  793.    float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
  794.    int incX = -1;
  795.    float x_expected[] = { 1.05233f, 0.79657f, -0.566883f, 1.46031f, -0.437f, 0.592f };
  796.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  797.    {
  798.      int i;
  799.      for (i = 0; i < 3; i++) {
  800.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1265) real");
  801.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1265) imag");
  802.      };
  803.    };
  804.   };
  805.  
  806.  
  807.   {
  808.    int order = 102;
  809.    int trans = 111;
  810.    int uplo = 121;
  811.    int diag = 131;
  812.    int N = 3;
  813.    int K = 1;
  814.    int lda = 3;
  815.    float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
  816.    float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
  817.    int incX = -1;
  818.    float x_expected[] = { -0.735844f, 1.11782f, -0.28244f, 1.16117f, -0.66707f, 0.938302f };
  819.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  820.    {
  821.      int i;
  822.      for (i = 0; i < 3; i++) {
  823.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1266) real");
  824.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1266) imag");
  825.      };
  826.    };
  827.   };
  828.  
  829.  
  830.   {
  831.    int order = 102;
  832.    int trans = 111;
  833.    int uplo = 121;
  834.    int diag = 132;
  835.    int N = 3;
  836.    int K = 1;
  837.    int lda = 3;
  838.    float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
  839.    float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
  840.    int incX = -1;
  841.    float x_expected[] = { 0.11f, 0.787f, -0.406239f, 0.580226f, -0.171935f, 1.2125f };
  842.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  843.    {
  844.      int i;
  845.      for (i = 0; i < 3; i++) {
  846.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1267) real");
  847.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1267) imag");
  848.      };
  849.    };
  850.   };
  851.  
  852.  
  853.   {
  854.    int order = 102;
  855.    int trans = 111;
  856.    int uplo = 122;
  857.    int diag = 131;
  858.    int N = 3;
  859.    int K = 1;
  860.    int lda = 3;
  861.    float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
  862.    float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
  863.    int incX = -1;
  864.    float x_expected[] = { 1.70081f, 2.20477f, 1.32753f, -0.522112f, 0.0223652f, -0.62248f };
  865.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  866.    {
  867.      int i;
  868.      for (i = 0; i < 3; i++) {
  869.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1268) real");
  870.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1268) imag");
  871.      };
  872.    };
  873.   };
  874.  
  875.  
  876.   {
  877.    int order = 102;
  878.    int trans = 111;
  879.    int uplo = 122;
  880.    int diag = 132;
  881.    int N = 3;
  882.    int K = 1;
  883.    int lda = 3;
  884.    float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
  885.    float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
  886.    int incX = -1;
  887.    float x_expected[] = { 0.967596f, 0.693563f, -1.04022f, -0.09269f, -0.437f, 0.592f };
  888.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  889.    {
  890.      int i;
  891.      for (i = 0; i < 3; i++) {
  892.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1269) real");
  893.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1269) imag");
  894.      };
  895.    };
  896.   };
  897.  
  898.  
  899.   {
  900.    int order = 101;
  901.    int trans = 112;
  902.    int uplo = 121;
  903.    int diag = 131;
  904.    int N = 3;
  905.    int K = 1;
  906.    int lda = 3;
  907.    float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
  908.    float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
  909.    int incX = -1;
  910.    float x_expected[] = { -1.11985f, 0.801655f, 0.273814f, -1.09438f, -0.52531f, 0.166748f };
  911.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  912.    {
  913.      int i;
  914.      for (i = 0; i < 3; i++) {
  915.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1270) real");
  916.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1270) imag");
  917.      };
  918.    };
  919.   };
  920.  
  921.  
  922.   {
  923.    int order = 101;
  924.    int trans = 112;
  925.    int uplo = 121;
  926.    int diag = 132;
  927.    int N = 3;
  928.    int K = 1;
  929.    int lda = 3;
  930.    float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
  931.    float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
  932.    int incX = -1;
  933.    float x_expected[] = { 0.266087f, 0.618557f, 0.031897f, -0.914419f, -0.134f, 0.179f };
  934.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  935.    {
  936.      int i;
  937.      for (i = 0; i < 3; i++) {
  938.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1271) real");
  939.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1271) imag");
  940.      };
  941.    };
  942.   };
  943.  
  944.  
  945.   {
  946.    int order = 101;
  947.    int trans = 112;
  948.    int uplo = 122;
  949.    int diag = 131;
  950.    int N = 3;
  951.    int K = 1;
  952.    int lda = 3;
  953.    float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
  954.    float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
  955.    int incX = -1;
  956.    float x_expected[] = { -0.762749f, -0.016292f, 1.59299f, 0.158751f, -4.75603f, -1.78591f };
  957.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  958.    {
  959.      int i;
  960.      for (i = 0; i < 3; i++) {
  961.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1272) real");
  962.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1272) imag");
  963.      };
  964.    };
  965.   };
  966.  
  967.  
  968.   {
  969.    int order = 101;
  970.    int trans = 112;
  971.    int uplo = 122;
  972.    int diag = 132;
  973.    int N = 3;
  974.    int K = 1;
  975.    int lda = 3;
  976.    float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
  977.    float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
  978.    int incX = -1;
  979.    float x_expected[] = { -0.509f, 0.608f, -0.332731f, -1.24444f, 0.262904f, 1.21961f };
  980.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  981.    {
  982.      int i;
  983.      for (i = 0; i < 3; i++) {
  984.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1273) real");
  985.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1273) imag");
  986.      };
  987.    };
  988.   };
  989.  
  990.  
  991.   {
  992.    int order = 102;
  993.    int trans = 112;
  994.    int uplo = 121;
  995.    int diag = 131;
  996.    int N = 3;
  997.    int K = 1;
  998.    int lda = 3;
  999.    float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
  1000.    float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
  1001.    int incX = -1;
  1002.    float x_expected[] = { -1.76046f, 0.0455463f, 1.38348f, 0.700097f, -0.669451f, 0.321896f };
  1003.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1004.    {
  1005.      int i;
  1006.      for (i = 0; i < 3; i++) {
  1007.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1274) real");
  1008.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1274) imag");
  1009.      };
  1010.    };
  1011.   };
  1012.  
  1013.  
  1014.   {
  1015.    int order = 102;
  1016.    int trans = 112;
  1017.    int uplo = 121;
  1018.    int diag = 132;
  1019.    int N = 3;
  1020.    int K = 1;
  1021.    int lda = 3;
  1022.    float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
  1023.    float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
  1024.    int incX = -1;
  1025.    float x_expected[] = { 0.151523f, 0.78611f, 0.120309f, -1.01387f, -0.134f, 0.179f };
  1026.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1027.    {
  1028.      int i;
  1029.      for (i = 0; i < 3; i++) {
  1030.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1275) real");
  1031.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1275) imag");
  1032.      };
  1033.    };
  1034.   };
  1035.  
  1036.  
  1037.   {
  1038.    int order = 102;
  1039.    int trans = 112;
  1040.    int uplo = 122;
  1041.    int diag = 131;
  1042.    int N = 3;
  1043.    int K = 1;
  1044.    int lda = 3;
  1045.    float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
  1046.    float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
  1047.    int incX = -1;
  1048.    float x_expected[] = { -1.00779f, -0.620278f, 0.81164f, -1.90759f, -1.32022f, 1.48356f };
  1049.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1050.    {
  1051.      int i;
  1052.      for (i = 0; i < 3; i++) {
  1053.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1276) real");
  1054.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1276) imag");
  1055.      };
  1056.    };
  1057.   };
  1058.  
  1059.  
  1060.   {
  1061.    int order = 102;
  1062.    int trans = 112;
  1063.    int uplo = 122;
  1064.    int diag = 132;
  1065.    int N = 3;
  1066.    int K = 1;
  1067.    int lda = 3;
  1068.    float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
  1069.    float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
  1070.    int incX = -1;
  1071.    float x_expected[] = { -0.509f, 0.608f, -0.503138f, -1.26818f, 0.176615f, 0.447668f };
  1072.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1073.    {
  1074.      int i;
  1075.      for (i = 0; i < 3; i++) {
  1076.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1277) real");
  1077.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1277) imag");
  1078.      };
  1079.    };
  1080.   };
  1081.  
  1082.  
  1083.   {
  1084.    int order = 101;
  1085.    int trans = 113;
  1086.    int uplo = 121;
  1087.    int diag = 131;
  1088.    int N = 3;
  1089.    int K = 1;
  1090.    int lda = 3;
  1091.    float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
  1092.    float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
  1093.    int incX = -1;
  1094.    float x_expected[] = { -0.613838f, -1.13321f, -1.34847f, 0.0432903f, 0.0879552f, -0.479334f };
  1095.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1096.    {
  1097.      int i;
  1098.      for (i = 0; i < 3; i++) {
  1099.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1278) real");
  1100.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1278) imag");
  1101.      };
  1102.    };
  1103.   };
  1104.  
  1105.  
  1106.   {
  1107.    int order = 101;
  1108.    int trans = 113;
  1109.    int uplo = 121;
  1110.    int diag = 132;
  1111.    int N = 3;
  1112.    int K = 1;
  1113.    int lda = 3;
  1114.    float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
  1115.    float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
  1116.    int incX = -1;
  1117.    float x_expected[] = { 0.76323f, -1.23595f, 0.943058f, -0.618694f, 0.296f, 0.034f };
  1118.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1119.    {
  1120.      int i;
  1121.      for (i = 0; i < 3; i++) {
  1122.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1279) real");
  1123.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1279) imag");
  1124.      };
  1125.    };
  1126.   };
  1127.  
  1128.  
  1129.   {
  1130.    int order = 101;
  1131.    int trans = 113;
  1132.    int uplo = 122;
  1133.    int diag = 131;
  1134.    int N = 3;
  1135.    int K = 1;
  1136.    int lda = 3;
  1137.    float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
  1138.    float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
  1139.    int incX = -1;
  1140.    float x_expected[] = { -1.15557f, -2.50103f, -3.85402f, -1.04833f, 0.414582f, 5.91218f };
  1141.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1142.    {
  1143.      int i;
  1144.      for (i = 0; i < 3; i++) {
  1145.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1280) real");
  1146.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1280) imag");
  1147.      };
  1148.    };
  1149.   };
  1150.  
  1151.  
  1152.   {
  1153.    int order = 101;
  1154.    int trans = 113;
  1155.    int uplo = 122;
  1156.    int diag = 132;
  1157.    int N = 3;
  1158.    int K = 1;
  1159.    int lda = 3;
  1160.    float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
  1161.    float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
  1162.    int incX = -1;
  1163.    float x_expected[] = { -0.037f, -0.599f, 1.39953f, -0.064424f, 1.0801f, -0.481747f };
  1164.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1165.    {
  1166.      int i;
  1167.      for (i = 0; i < 3; i++) {
  1168.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1281) real");
  1169.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1281) imag");
  1170.      };
  1171.    };
  1172.   };
  1173.  
  1174.  
  1175.   {
  1176.    int order = 102;
  1177.    int trans = 113;
  1178.    int uplo = 121;
  1179.    int diag = 131;
  1180.    int N = 3;
  1181.    int K = 1;
  1182.    int lda = 3;
  1183.    float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
  1184.    float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
  1185.    int incX = -1;
  1186.    float x_expected[] = { -3.0802f, -9.09377f, -1.05845f, 0.99239f, 0.259763f, -0.687744f };
  1187.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1188.    {
  1189.      int i;
  1190.      for (i = 0; i < 3; i++) {
  1191.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1282) real");
  1192.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1282) imag");
  1193.      };
  1194.    };
  1195.   };
  1196.  
  1197.  
  1198.   {
  1199.    int order = 102;
  1200.    int trans = 113;
  1201.    int uplo = 121;
  1202.    int diag = 132;
  1203.    int N = 3;
  1204.    int K = 1;
  1205.    int lda = 3;
  1206.    float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
  1207.    float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
  1208.    int incX = -1;
  1209.    float x_expected[] = { -0.513897f, 0.632031f, 1.14112f, -0.580648f, 0.296f, 0.034f };
  1210.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1211.    {
  1212.      int i;
  1213.      for (i = 0; i < 3; i++) {
  1214.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1283) real");
  1215.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1283) imag");
  1216.      };
  1217.    };
  1218.   };
  1219.  
  1220.  
  1221.   {
  1222.    int order = 102;
  1223.    int trans = 113;
  1224.    int uplo = 122;
  1225.    int diag = 131;
  1226.    int N = 3;
  1227.    int K = 1;
  1228.    int lda = 3;
  1229.    float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
  1230.    float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
  1231.    int incX = -1;
  1232.    float x_expected[] = { 0.360899f, -0.456643f, -2.31803f, 0.257877f, 1.56928f, -0.922115f };
  1233.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1234.    {
  1235.      int i;
  1236.      for (i = 0; i < 3; i++) {
  1237.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1284) real");
  1238.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1284) imag");
  1239.      };
  1240.    };
  1241.   };
  1242.  
  1243.  
  1244.   {
  1245.    int order = 102;
  1246.    int trans = 113;
  1247.    int uplo = 122;
  1248.    int diag = 132;
  1249.    int N = 3;
  1250.    int K = 1;
  1251.    int lda = 3;
  1252.    float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
  1253.    float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
  1254.    int incX = -1;
  1255.    float x_expected[] = { -0.037f, -0.599f, 0.875872f, -1.03683f, -0.198184f, -0.207572f };
  1256.    cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1257.    {
  1258.      int i;
  1259.      for (i = 0; i < 3; i++) {
  1260.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1285) real");
  1261.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1285) imag");
  1262.      };
  1263.    };
  1264.   };
  1265.  
  1266.  
  1267.   {
  1268.    int order = 101;
  1269.    int trans = 111;
  1270.    int uplo = 121;
  1271.    int diag = 131;
  1272.    int N = 3;
  1273.    int K = 1;
  1274.    int lda = 3;
  1275.    double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
  1276.    double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
  1277.    int incX = -1;
  1278.    double x_expected[] = { 0.0490338308139, -0.158433417494, 0.261604043488, 1.28058846321, 1.77633350191, -1.07039599422 };
  1279.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1280.    {
  1281.      int i;
  1282.      for (i = 0; i < 3; i++) {
  1283.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1286) real");
  1284.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1286) imag");
  1285.      };
  1286.    };
  1287.   };
  1288.  
  1289.  
  1290.   {
  1291.    int order = 101;
  1292.    int trans = 111;
  1293.    int uplo = 121;
  1294.    int diag = 132;
  1295.    int N = 3;
  1296.    int K = 1;
  1297.    int lda = 3;
  1298.    double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
  1299.    double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
  1300.    int incX = -1;
  1301.    double x_expected[] = { -0.123, 0.122, 0.96534, 0.346049, 1.067212328, 0.445330131 };
  1302.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1303.    {
  1304.      int i;
  1305.      for (i = 0; i < 3; i++) {
  1306.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1287) real");
  1307.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1287) imag");
  1308.      };
  1309.    };
  1310.   };
  1311.  
  1312.  
  1313.   {
  1314.    int order = 101;
  1315.    int trans = 111;
  1316.    int uplo = 122;
  1317.    int diag = 131;
  1318.    int N = 3;
  1319.    int K = 1;
  1320.    int lda = 3;
  1321.    double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
  1322.    double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
  1323.    int incX = -1;
  1324.    double x_expected[] = { 72.7437666278, 10.4206532927, -4.34946941374, -14.8012581742, 2.01859491883, -1.53922125931 };
  1325.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1326.    {
  1327.      int i;
  1328.      for (i = 0; i < 3; i++) {
  1329.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1288) real");
  1330.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1288) imag");
  1331.      };
  1332.    };
  1333.   };
  1334.  
  1335.  
  1336.   {
  1337.    int order = 101;
  1338.    int trans = 111;
  1339.    int uplo = 122;
  1340.    int diag = 132;
  1341.    int N = 3;
  1342.    int K = 1;
  1343.    int lda = 3;
  1344.    double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
  1345.    double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
  1346.    int incX = -1;
  1347.    double x_expected[] = { -0.464775024, 0.662224708, -0.0457, 0.610264, 0.942, 0.98 };
  1348.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1349.    {
  1350.      int i;
  1351.      for (i = 0; i < 3; i++) {
  1352.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1289) real");
  1353.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1289) imag");
  1354.      };
  1355.    };
  1356.   };
  1357.  
  1358.  
  1359.   {
  1360.    int order = 102;
  1361.    int trans = 111;
  1362.    int uplo = 121;
  1363.    int diag = 131;
  1364.    int N = 3;
  1365.    int K = 1;
  1366.    int lda = 3;
  1367.    double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
  1368.    double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
  1369.    int incX = -1;
  1370.    double x_expected[] = { -0.591747295323, -0.534096923761, -4.60251824353, 1.70172936273, -4.94687072873, -3.32536493524 };
  1371.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1372.    {
  1373.      int i;
  1374.      for (i = 0; i < 3; i++) {
  1375.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1290) real");
  1376.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1290) imag");
  1377.      };
  1378.    };
  1379.   };
  1380.  
  1381.  
  1382.   {
  1383.    int order = 102;
  1384.    int trans = 111;
  1385.    int uplo = 121;
  1386.    int diag = 132;
  1387.    int N = 3;
  1388.    int K = 1;
  1389.    int lda = 3;
  1390.    double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
  1391.    double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
  1392.    int incX = -1;
  1393.    double x_expected[] = { -0.123, 0.122, 0.807692, 0.373091, 0.384974988, 1.400879194 };
  1394.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1395.    {
  1396.      int i;
  1397.      for (i = 0; i < 3; i++) {
  1398.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1291) real");
  1399.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1291) imag");
  1400.      };
  1401.    };
  1402.   };
  1403.  
  1404.  
  1405.   {
  1406.    int order = 102;
  1407.    int trans = 111;
  1408.    int uplo = 122;
  1409.    int diag = 131;
  1410.    int N = 3;
  1411.    int K = 1;
  1412.    int lda = 3;
  1413.    double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
  1414.    double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
  1415.    int incX = -1;
  1416.    double x_expected[] = { -0.129998778267, -0.116630230861, 0.993340886904, 0.530739563688, 1.55891621291, -0.284019181928 };
  1417.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1418.    {
  1419.      int i;
  1420.      for (i = 0; i < 3; i++) {
  1421.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1292) real");
  1422.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1292) imag");
  1423.      };
  1424.    };
  1425.   };
  1426.  
  1427.  
  1428.   {
  1429.    int order = 102;
  1430.    int trans = 111;
  1431.    int uplo = 122;
  1432.    int diag = 132;
  1433.    int N = 3;
  1434.    int K = 1;
  1435.    int lda = 3;
  1436.    double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
  1437.    double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
  1438.    int incX = -1;
  1439.    double x_expected[] = { 0.107496032, 0.025821594, 1.444898, -0.239924, 0.942, 0.98 };
  1440.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1441.    {
  1442.      int i;
  1443.      for (i = 0; i < 3; i++) {
  1444.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1293) real");
  1445.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1293) imag");
  1446.      };
  1447.    };
  1448.   };
  1449.  
  1450.  
  1451.   {
  1452.    int order = 101;
  1453.    int trans = 112;
  1454.    int uplo = 121;
  1455.    int diag = 131;
  1456.    int N = 3;
  1457.    int K = 1;
  1458.    int lda = 3;
  1459.    double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
  1460.    double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
  1461.    int incX = -1;
  1462.    double x_expected[] = { -0.825842176606, 0.212941473892, -0.548817434511, -0.703261551538, 0.0746069436827, 0.425751789407 };
  1463.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1464.    {
  1465.      int i;
  1466.      for (i = 0; i < 3; i++) {
  1467.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1294) real");
  1468.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1294) imag");
  1469.      };
  1470.    };
  1471.   };
  1472.  
  1473.  
  1474.   {
  1475.    int order = 101;
  1476.    int trans = 112;
  1477.    int uplo = 121;
  1478.    int diag = 132;
  1479.    int N = 3;
  1480.    int K = 1;
  1481.    int lda = 3;
  1482.    double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
  1483.    double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
  1484.    int incX = -1;
  1485.    double x_expected[] = { -0.619710352, 0.018225936, 1.211252, 0.891864, 0.293, -0.434 };
  1486.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1487.    {
  1488.      int i;
  1489.      for (i = 0; i < 3; i++) {
  1490.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1295) real");
  1491.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1295) imag");
  1492.      };
  1493.    };
  1494.   };
  1495.  
  1496.  
  1497.   {
  1498.    int order = 101;
  1499.    int trans = 112;
  1500.    int uplo = 122;
  1501.    int diag = 131;
  1502.    int N = 3;
  1503.    int K = 1;
  1504.    int lda = 3;
  1505.    double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
  1506.    double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
  1507.    int incX = -1;
  1508.    double x_expected[] = { 0.203289119964, 1.58288482537, -1.7720160159, 0.479463518178, -0.511241930019, -1.79333888299 };
  1509.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1510.    {
  1511.      int i;
  1512.      for (i = 0; i < 3; i++) {
  1513.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1296) real");
  1514.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1296) imag");
  1515.      };
  1516.    };
  1517.   };
  1518.  
  1519.  
  1520.   {
  1521.    int order = 101;
  1522.    int trans = 112;
  1523.    int uplo = 122;
  1524.    int diag = 132;
  1525.    int N = 3;
  1526.    int K = 1;
  1527.    int lda = 3;
  1528.    double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
  1529.    double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
  1530.    int incX = -1;
  1531.    double x_expected[] = { -0.373, 0.566, 0.618602, -0.084689, 0.887531803, -0.570220771 };
  1532.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1533.    {
  1534.      int i;
  1535.      for (i = 0; i < 3; i++) {
  1536.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1297) real");
  1537.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1297) imag");
  1538.      };
  1539.    };
  1540.   };
  1541.  
  1542.  
  1543.   {
  1544.    int order = 102;
  1545.    int trans = 112;
  1546.    int uplo = 121;
  1547.    int diag = 131;
  1548.    int N = 3;
  1549.    int K = 1;
  1550.    int lda = 3;
  1551.    double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
  1552.    double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
  1553.    int incX = -1;
  1554.    double x_expected[] = { 1.72799012007, 13.4612400765, 4.46126528205, -0.0212528722047, 0.627282377919, 0.302760084926 };
  1555.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1556.    {
  1557.      int i;
  1558.      for (i = 0; i < 3; i++) {
  1559.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1298) real");
  1560.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1298) imag");
  1561.      };
  1562.    };
  1563.   };
  1564.  
  1565.  
  1566.   {
  1567.    int order = 102;
  1568.    int trans = 112;
  1569.    int uplo = 121;
  1570.    int diag = 132;
  1571.    int N = 3;
  1572.    int K = 1;
  1573.    int lda = 3;
  1574.    double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
  1575.    double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
  1576.    int incX = -1;
  1577.    double x_expected[] = { -1.280839615, 1.560525655, 1.167331, 0.179227, 0.293, -0.434 };
  1578.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1579.    {
  1580.      int i;
  1581.      for (i = 0; i < 3; i++) {
  1582.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1299) real");
  1583.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1299) imag");
  1584.      };
  1585.    };
  1586.   };
  1587.  
  1588.  
  1589.   {
  1590.    int order = 102;
  1591.    int trans = 112;
  1592.    int uplo = 122;
  1593.    int diag = 131;
  1594.    int N = 3;
  1595.    int K = 1;
  1596.    int lda = 3;
  1597.    double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
  1598.    double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
  1599.    int incX = -1;
  1600.    double x_expected[] = { -0.594503951847, 0.00287302167266, -1.08185265666, -0.859860374254, 0.0331027077244, 1.28233265933 };
  1601.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1602.    {
  1603.      int i;
  1604.      for (i = 0; i < 3; i++) {
  1605.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1300) real");
  1606.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1300) imag");
  1607.      };
  1608.    };
  1609.   };
  1610.  
  1611.  
  1612.   {
  1613.    int order = 102;
  1614.    int trans = 112;
  1615.    int uplo = 122;
  1616.    int diag = 132;
  1617.    int N = 3;
  1618.    int K = 1;
  1619.    int lda = 3;
  1620.    double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
  1621.    double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
  1622.    int incX = -1;
  1623.    double x_expected[] = { -0.373, 0.566, 1.16074, 0.50314, -0.20669608, 0.37525144 };
  1624.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1625.    {
  1626.      int i;
  1627.      for (i = 0; i < 3; i++) {
  1628.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1301) real");
  1629.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1301) imag");
  1630.      };
  1631.    };
  1632.   };
  1633.  
  1634.  
  1635.   {
  1636.    int order = 101;
  1637.    int trans = 113;
  1638.    int uplo = 121;
  1639.    int diag = 131;
  1640.    int N = 3;
  1641.    int K = 1;
  1642.    int lda = 3;
  1643.    double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
  1644.    double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
  1645.    int incX = -1;
  1646.    double x_expected[] = { 0.0654496252357, 0.224007771015, -0.752486084395, -0.554870892947, -0.587163401057, 0.166737652215 };
  1647.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1648.    {
  1649.      int i;
  1650.      for (i = 0; i < 3; i++) {
  1651.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1302) real");
  1652.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1302) imag");
  1653.      };
  1654.    };
  1655.   };
  1656.  
  1657.  
  1658.   {
  1659.    int order = 101;
  1660.    int trans = 113;
  1661.    int uplo = 121;
  1662.    int diag = 132;
  1663.    int N = 3;
  1664.    int K = 1;
  1665.    int lda = 3;
  1666.    double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
  1667.    double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
  1668.    int incX = -1;
  1669.    double x_expected[] = { -0.595558802, -1.147174647, 0.589506, -0.500919, -0.126, 0.459 };
  1670.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1671.    {
  1672.      int i;
  1673.      for (i = 0; i < 3; i++) {
  1674.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1303) real");
  1675.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1303) imag");
  1676.      };
  1677.    };
  1678.   };
  1679.  
  1680.  
  1681.   {
  1682.    int order = 101;
  1683.    int trans = 113;
  1684.    int uplo = 122;
  1685.    int diag = 131;
  1686.    int N = 3;
  1687.    int K = 1;
  1688.    int lda = 3;
  1689.    double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
  1690.    double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
  1691.    int incX = -1;
  1692.    double x_expected[] = { 3.39346077201, 0.652889512141, -2.33602680355, -2.7859245153, -5.04672104102, -0.334110541026 };
  1693.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1694.    {
  1695.      int i;
  1696.      for (i = 0; i < 3; i++) {
  1697.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1304) real");
  1698.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1304) imag");
  1699.      };
  1700.    };
  1701.   };
  1702.  
  1703.  
  1704.   {
  1705.    int order = 101;
  1706.    int trans = 113;
  1707.    int uplo = 122;
  1708.    int diag = 132;
  1709.    int N = 3;
  1710.    int K = 1;
  1711.    int lda = 3;
  1712.    double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
  1713.    double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
  1714.    int incX = -1;
  1715.    double x_expected[] = { 0.028, -0.804, -0.109456, -0.217192, -0.41110804, 0.41693792 };
  1716.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1717.    {
  1718.      int i;
  1719.      for (i = 0; i < 3; i++) {
  1720.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1305) real");
  1721.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1305) imag");
  1722.      };
  1723.    };
  1724.   };
  1725.  
  1726.  
  1727.   {
  1728.    int order = 102;
  1729.    int trans = 113;
  1730.    int uplo = 121;
  1731.    int diag = 131;
  1732.    int N = 3;
  1733.    int K = 1;
  1734.    int lda = 3;
  1735.    double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
  1736.    double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
  1737.    int incX = -1;
  1738.    double x_expected[] = { 7.16970224467, -0.772071373678, 0.833386981173, -0.673826630129, -0.26524050899, 0.465327628365 };
  1739.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1740.    {
  1741.      int i;
  1742.      for (i = 0; i < 3; i++) {
  1743.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1306) real");
  1744.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1306) imag");
  1745.      };
  1746.    };
  1747.   };
  1748.  
  1749.  
  1750.   {
  1751.    int order = 102;
  1752.    int trans = 113;
  1753.    int uplo = 121;
  1754.    int diag = 132;
  1755.    int N = 3;
  1756.    int K = 1;
  1757.    int lda = 3;
  1758.    double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
  1759.    double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
  1760.    int incX = -1;
  1761.    double x_expected[] = { 0.471459157, -1.566755859, 0.940839, 0.357132, -0.126, 0.459 };
  1762.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1763.    {
  1764.      int i;
  1765.      for (i = 0; i < 3; i++) {
  1766.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1307) real");
  1767.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1307) imag");
  1768.      };
  1769.    };
  1770.   };
  1771.  
  1772.  
  1773.   {
  1774.    int order = 102;
  1775.    int trans = 113;
  1776.    int uplo = 122;
  1777.    int diag = 131;
  1778.    int N = 3;
  1779.    int K = 1;
  1780.    int lda = 3;
  1781.    double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
  1782.    double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
  1783.    int incX = -1;
  1784.    double x_expected[] = { -0.909961830373, 0.118063054039, -0.0169425582229, -1.00055409731, -1.37205489923, 0.994032418785 };
  1785.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1786.    {
  1787.      int i;
  1788.      for (i = 0; i < 3; i++) {
  1789.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1308) real");
  1790.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1308) imag");
  1791.      };
  1792.    };
  1793.   };
  1794.  
  1795.  
  1796.   {
  1797.    int order = 102;
  1798.    int trans = 113;
  1799.    int uplo = 122;
  1800.    int diag = 132;
  1801.    int N = 3;
  1802.    int K = 1;
  1803.    int lda = 3;
  1804.    double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
  1805.    double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
  1806.    int incX = -1;
  1807.    double x_expected[] = { 0.028, -0.804, -0.118596, 0.160828, -0.059271004, 0.294435972 };
  1808.    cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1809.    {
  1810.      int i;
  1811.      for (i = 0; i < 3; i++) {
  1812.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1309) real");
  1813.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1309) imag");
  1814.      };
  1815.    };
  1816.   };
  1817.  
  1818.  
  1819. }
  1820.